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ABSTRACT 

We investigate the predicted present-day temperature profiles of the hot, X-ray emit- 
ting gas in galaxy clusters for two cosmological models - a current best-guess ACDM 
model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" 
of clusters are derived from high-resolution (15 h~^kpc) simulations which make use of a 
sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully 
captures the shocks which are essential for correctly modelling cluster temperatures. 
We show that the temperature structure on Mpc-scales is highly complex and non- 
isothermal. However, the temperature profiles of the simulated ACDM and SCDM 
clusters are remarkably similar and drop-off as T oc (1 + r/a^)"^ where ~ ryir/l.h 
and 5 ~ 1.6. This decrease is in good agreement with the observational results of 
Markevitch et al. (1998) but diverges, primarily in the innermost regions, from their 
fit which assumes a polytropic equation of state. Our result is also in good agreement 
with a recent sample of clusters observed by BeppoSKY^ though there is some indication 
of missing physics at small radii (r < 0.2 Vyir). We discuss the interpretation of our 
results and make predictions for new x-ray observations that will extend to larger radii 
than previously possible. Finally, we show that, for r > 0.2 r^j^, our universal tempera- 
ture profile is consistent with our most recent simulations which include both radiative 
cooling and supernovae feedback. 

Subject headings: cosmology: theory - galaxies: clusters: general - hydrodynamics - intergalactic 
medium - X-rays:galaxies 



1. Introduction 

X-ray based cluster mass determinations 
generally assume that morphologically-symmetric 
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(i.e. non-merging) clusters are isothermal. 
However, this accepted wisdom has recently 
been challenged by Markevitch et al. (1998; 
hereafter M98) who found evidence for de- 
creasing temperature profiles in a sample of 
nearby hot clusters (> 3.5 keV) observed with 
ASCA. A subsample of 17 regular /symmetric 
clusters displayed remarkably similar temper- 
ature profiles (when normalized and scaled 
by the virial radius) consistent with T oc 
[1 + (r/rc)2]-3/3(T-i)/2 where 7 = 1.24^:20 
and (3 = 2/3. The typical decrease is there- 
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fore a factor of ~2 in going from 1 to 6 core 
radii (or .09 to 0.5 virial radii). This result re- 
mains controversial as three subsequent stud- 
ies of large samples of clusters concluded that 
the majority of cluster temperature profiles 
show little, or no, decrease with radius (Irwin, 
Bregman, & Evrard 1999; White 2000; Irwin 
&; Bregman 2000). Most recently, De Grandi 
&; Molendi (2002) have presented a composite 
temperature profile based on BeppoSAji data 
which exhibits an isothermal core and then 
decreases quickly. Here we present results de- 
rived from recent high-resolution numerical 
simulations and show that there appears to 
be a universal temperature profile which de- 
clines significantly even within half a virial 
radius of the cluster center. We also show 
that these simulated profiles arc consistent 
with the most recent cluster observations ex- 
cept in the very innermost regions of a cluster 
(r < 0.2 ryir). 

2. Simulations 

The simulated clusters used in this work 
are part of an ongoing project to create ex- 
tensive numerical "catalogs" of clusters. This 
catalog has already been used to investigate 
the effect of large-scale structure on weak 
gravitational lensing mass estimates (Metzler 
et al. 1999, 2001) and the role of "temper- 
ature bias" in Sunyaev-Zeldovich determina- 
tions of the Hubble constant (Lin, Norman, 
& Bryan 2001). The cluster simulation data 
are available on-line at the Simulated Cluster 
Archive: http://sca.ncsa.uiuc.edu (Norman et 
al. 1999). 

The cosmological model adopted for the 
bulk of the simulations is a current "best- 
guess" fiat universe with a cosmological 
constant (ACDM) and parameters: f^o = 
0.3, = 0.026, nA = 0.7, h = 0.7 and 
as = 0.928. The 25 most-massive clusters in 
the 256 h~^ Mpc volume have been fully ana- 
lyzed and comprise the ACDM sample used in 



this work. For comparison, a smaller sample 
of 10 SCDM clusters has been simulated in the 
same volume (and with the same random-seed 
for the initial fluctuations) with parameters: 
Oo = 1, Ob = 0.076, , ^ = 0.5 and as = 0.6 

An initial low-resolution simulation was 
used to identify clusters in the 256 h~^ Mpc 
volume and then high-resolution simulations 
were performed which evolved the entire vol- 
ume but adaptively refined the region from 
which a specific cluster formed. The root, or 
level (LO) , grid covers the entire periodic do- 
main with 128^ cells. The Lagrangian volume 
of a given cluster is further refined with two 
levels of subgrids (LI and L2). each with twice 
the spatial resolution and eight times the mass 
resolution of the previous one, in order to re- 
fine the initial conditions within this volume. 
An additional five levels of refinement are in- 
troduced adaptively and automatically within 
the L2 grid as the simulation progresses. The 
clusters form from particles on the L2 grid 
which have a mass of 9 X lO^h-^M0. The 
best spatial resolution is 15.6 h~^kpc. 

The cosmological adaptive-mesh refine- 
ment (AMR) code used here is described 
in detail in Bryan (1999) and Norman & 
Bryan (1999). Briefly, this code solves the 
equations of hydrodynamics using the higher- 
order accurate piecewise parabolic method 
(PPM; Collela & Woodward 1984) and em- 
ploys an adaptive particle-mesh algorithm 
(using second-order accurate TSC interpo- 
lation) for the dark matter particles. The 
Eulerian hydrodynamics algorithm exhibits 
several crucial characteristics including excel- 
lent shock-capturing (within 1-2 zones with 
correct entropy generation), and equal accu- 
racy in high or low density regions (for the 
same zone size). The key component of the 
code is its use of structured adaptive mesh re- 
finement (Berger Sz Collela 1989) in order to 
enhance resolution when and where needed. 

Our code was tested as part of the Santa 
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Barbara cluster comparison project (Frenk et 
al. 1999) in which 12 groups simulated a 
Coma-like cluster using a variety of codes 
and resolutions Our results (curves labeled 
"bryan") are among the highest resolution re- 
sults presented in that paper (central reso- 
lution of 7.8h~^ kpc). We agree well with 
the highest resolution AP^MSPH results ex- 
cept that our central temperature is some- 
what higher, and our radial temperature pro- 
files continue to rise at small radii rather than 
flattening within 200 kpc. Frenk et al. (1999) 
note that this discrepancy is more or less 
present in all the grid-based codes but pro- 
vide no explanation. To check whether this 
might be a result of numerical resolution or al- 
gorithms employed, wc recomputed the Santa 
Barbara cluster with low, medium and high 
resolution initial conditions corresponding to 
64^,1283 and 256^ particles/cells in the LO 
grid, and L7 resolutions equal to 15.6. 7.8, 
and 3.9h^^ kpc, respectively. Instead of PPM, 
we also tested the artificial viscosity, finite- 
difference algorithm for gas dynamics used 
in the ZEUS codes (Stone & Norman 1992). 
In all cases we found the temperature profile 
agreed with our result presented in Frenk et 
al. to within a few percent at all radii; i.e., 
our temperature profile has converged. Note 
also that there is excellent agreement between 
our results for the Santa Barbara cluster and 
those of a new, completely independent AMR 
code (see Fig.l of Kravtsov, Klypin, & Hoff- 
man 2001). 

3. Results 

For each simulation, clusters were identi- 
fied in the refined subvolume and their virial 
radii (ivir) determined as the radius at which 
the mean enclosed density was Ac times the 
critical density (at z = 0, Ac ~ 101 for 
ACDM and 178 for SCDM; e.g., Bryan & 
Norman 1998). Projected gas density, X-ray 
surface brightness, and emission-weighted gas 
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Fig. 1. — Histograms of the temperature, 
mass and virial radius distributions for the 
ACDM sample of simulated clusters. 

temperature images (5 h~^ Mpc on a side) 
were constructed for each cluster assuming 
a Raymond-Smith spectrum with metallicity 
0.3 solar in the (1.5-11.0) keV bandpass used 
in the analysis of M98. 2D radial temper- 
ature profiles (centered on the X-ray peak) 
were determined from the temperature im- 
ages and a global temperature (T^:) calculated 
by summing the emission-weighted tempera- 
ture within 1.0 h~^Mpc. This procedure was 
intended to duplicate, as closely as possible, 
that of M98. 

The distribution of cluster temperatures 
{Tx)i virial masses (M^j^), and virial radii 
(ryir) in the ACDM sample are shown in 
Fig. 1. The complex, non- isothermal struc- 
ture of the clusters is shown in the thumbnail 
images in Fig. 2. Although complicated, a 
number of patterns can be discerned: the core 
of the cluster is generally hotter than the out- 
skirts; low-temperature "holes" correspond to 
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Fig. 2. — Projected emission- weighted tem- 
perature maps of the ACDM clusters at 2; = 0. 
Each image is 5 h~^ Mpc on a side. Clusters 
are ordered from most massive (top left) to 
least massive (bottom right). 



inf ailing groups with lower masses and hence, 
lower temperatures. In addition, a few clus- 
ters show strong shocks due to recent major 
mergers, but these are generally not spher- 
ically symmetric and tend to be of order a 
Mpc from the core and so might be difficult 
to detect observationally. The evolution of the 
temperature structure and the significant dif- 
ferences between SCDM and ACDM clusters 
at z = 0.5 will be discussed in a forthcoming 
paper (Loken et al. 2002). 

The key result of this current paper con- 
cerns the 2D temperature profiles shown in 
Figs. 3 and 4. A given cluster's temperature 
profile is normalized by its global temperature 
(Tx) and the radius is scaled by the cluster's 
virial radius. Note that the resulting profiles 
are remarkably similar even when comparing 



ACDM and SCDM clusters. In addition, the 
majority of the profiles fall within the band 
drawn by M98 (reproduced in Figs. 3 and 4) 
to encompass the majority of their data points 
and error bars. 

Fig. 4 is a "tidied" version of Fig. 3 which 
includes only those clusters whose X-ray ap- 
pearance is fairly regular and symmetric. This 
is exactly the same selection cut imposed by 
M98 and others searching for a general tem- 
perature profile. However, even the excluded 
profiles are, for the most part, not dramati- 
cally different in general appearance. The dis- 
carded profiles generally have strong bumps 
or wiggles at the radius corresponding to a 
merging subcluster. 13 of 22 ACDM cluster 
profiles are kept in Fig. 4 (similarly, M98 find 
that 17 of their 33 clusters are symmetric). 
However, only 3 of 10 SCDM cluster are suf- 
ficiently regular; presumably this reflects the 
fact that clusters are still forming at 2: = in 
an SCDM universe. 

We find that the temperature profiles in 
Fig. 4 can be well fit with an expression of the 
form: T = To[l + r / ax]~^ , where Tq = 1.33, 
Ox = and 6 = 1.6 on the radial range 

(0.04-1.0) rmr- M98 chose instead to fit their 
temperature profiles with a polytropic equa- 
tion of state: T = To[l + (r/a^)2]-3/5(-T-i)/2. 
Using this form on the same radial range 
as M98, (0.09 - 0.55) r^ir, we find To = 
1.15, ax = rmr/7.2, and 7 = 1.26 (assum- 
ing /? = 2/3). M98 found a very similar fit: 
7 = 1.241 s-i^d Qx = r^ir /II- We note, how- 
ever, that the shape of the fit changes when 
the fitted range is extended to (0.04 — 1.0) ryir 
in which case: To = 1.17, ax = '^mr/5 and 
7 = 1.39. This sensitivity to the choice of 
fitting range appears to be due to the fact 
that the scatter in the data decreases for 
r >~ 0.6 r^ir and suggests that observations 
must extend to near the virial radius in order 
to converge on the true, universal temperature 
profile. 
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Fig. 3. — Temperature profiles for all the clus- 
ters in our ACDM (top) and SCDM (bottom) 
samples. Temperatures and radii are scaled in 
terms of the cluster's global temperature and 
virial radius, respectively. The heavy, black 
lines are reproduced from Fig. 7 of M98 and 
enclose most of their data points and error 
bars. 

Though the two forms for the temperature 
profile both yield acceptable fits, we prefer the 
first form for two reasons. First, the key dif- 
ference between the two fits is that the poly- 
tropic form turns over and begins to level-off 
at small radii while the first one continues to 
rise (just as in Fig. 4). If the fits are extended 
within r = 0.09 rmr, then the polytropic 
form yields an increasingly poor fit. It is dif- 
ficult to draw any strong conclusions, how- 
ever, because this is exactly the region within 
which non-adiabatic effects would become im- 
portant and, as well, our limiting resolution is 
~ 0.01 rmr for a typical cluster. 

The second argument against the poly- 
tropic form is the fact that it is not well- 
motivated physically. In fact, it can be de- 
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Fig. 4. — Temperature profiles for the sym- 
metric clusters in our ACDM (blue) and 
SCDM (red) samples. Scaling and the obser- 
vational bounds from M98 are as in Fig. 3. 

rived by assuming that the gas is: isother- 
mal, follows a /3-model for the X-ray surface 
brightness, and obeys the polytropic equation 
of state, T = ToP^~^ ■ Clearly, this is an in- 
consistent derivation. Moreover, the assump- 
tion of a polytropic equation of state amounts 
to assuming that the gas is isentropic which 
clearly contradicts the outwardly rising en- 
tropy profiles which we and others find from 
simulations (e.g., Frenk et al. 1999). 

Recently, De Grandi & Molendi (2002) 
have analyzed 5eppoSAX data for a sample 
of 21 clusters with and without cooling flows. 
They concluded that non-cooling flow cluster 
temperature profiles are flat in the innermost 
regions (r < 0.2 ryir) and that the profiles of 
both types of clusters decline rapidly with ra- 
dius outside r ~ 0.2 Vmr- Their mean data 
points for the two types of clusters are shown 
in Fig 5 along with profiles for two of our clus- 
ters which are discussed in §3.1. Again, the 
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Fig. 5. — Mean temperature profiles (with 
1-a errors) for a sample of 11 cooling flow 
(filled triangles) and 10 non-cooling flow clus- 
ters (filled circles) observed by BeppoSAK (De 
Grandi &; Molcndi 2002) along with our nu- 
merical results for a particular cluster simu- 
lated with and without radiative cooling. The 
grey bands represent the M98 observational 
bounds as in Fig. 3. 

shape and normalization of the observational 
profiles agree well with the simulations though 
the flat innermost region may signal the need 
for some additional physics in the simulations 
(e.g. De Grandi & Molendi 2002). 

3.1. Role of Additional Physics 

The present adiabatic simulations are the 
first step in a planned series of simulations 
which will incorporate non-adiabatic physical 
effects. Thus they form a baseline sample 
for evaluating the significance of additional 
physics. Though potentially important pro- 
cesses like cooling and galaxy feedback are not 
included in the simulations for the adiabatic 
cluster catalog, we would expect any effects 
to be strongly limited to the cluster core and 




0.2 



0.2 0.4 0.6 0.8 1.0 1.2 

Fig. 6. — Temperature profiles for a cluster 
from our adiabatic sample (solid line) and the 
same cluster rerun to include radiative cooling 
(dotted line). Our fit to the entire sample of 
symmetric clusters is shown as the dash-dot 
line. Also shown is the "Santa Barbara" clus- 
ter from a simulation incorporating cooling 
and supernovae feedback. See text for more 
details. The grey bands are again reproduced 
from Fig. 7 of M98 and enclose most of their 
data points and error bars. 

to have very little effect in the outer regions. 
These expectations now appear justified by 
our most recent non-adiabatic simulations. 

Two of the clusters in our sample have been 
re-run with identical parameters and resolu- 
tion in order to include the effects of radiative 
cooling (Motl ct al. 2002). These new runs are 
compared with the adiabatic simulations in 
Fig. 6. The temperature profiles of the cooling 
clusters are scaled as in the adiabatic case ex- 
cept that we omit the inner core (r < 0.1 r^r) 
when calculating the global temperature. In 
addition, one of us (MLN) has recently per- 
formed a higher resolution (7.8/i~^kpc) sim- 
ulation including both radiative cooling, and 
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Type II supernovae feedback. This cluster, 
which has the same initial parameters as the 
Santa Barbara comparison cluster (Frenk et 
al. 1999) is also included on the figure for 
comparison. Wc caution that this is just a 
first simulation and hence the full parameter 
space of star formation and feedback parame- 
ters has not been explored. 

Note that, as expected, the temperature 
profile of the cooling cluster is virtually indis- 
tinguishable from that of the adiabatic clus- 
ter outside of the cool core (r > 0.1 r^jr)- 
This strongly suggests that our conclusions 
about declining cluster temperature profiles 
are not the result of neglecting radiative cool- 
ing. Perhaps even more interesting is the 
higher-resolution simulation which incorpo- 
rates both cooling and feedback. Again, the 
shape of the corresponding temperature pro- 
file is remarkably similar to the other cases. 
Fig. 6 makes clear that the non-adiabatic 
profiles also fit well within the observational 
bounds determined by M98 

Electron thermal conduction, which is not 
included in our simulations, would have to be 
a significant fraction of its Spitzer value in 
order to erase temperature gradients on the 
scale we observe. The suppression of thermal 
conduction by tangled cluster magnetic fields 
is an open area of research which has, as yet, 
reached no firm conclusion (Chandran & Cow- 
ley, 1998; Malyshkin & Kulsrud, 2001). 

4. Discussion 

We have used our volume-limited catalogs 
of numerically-simulated clusters to demon- 
strate that both ACDM and SCDM clusters 
follow a universal cluster temperature pro- 
file with T/To = 1.3[1 + 1.5r/w]-i-6. This 
temperature profile agrees very well with the 
observationally-determined profile of M98 and 
the more recent i?eppoSAX data of DeGrandi 
& Molendi (2002) (Fig. 5). Our simulations 



also reveal a wealth of extremely complex and 
non-isothermal temperature structure (Fig. 2) 
which current X-ray telescopes may be able to 
probe. Interestingly, our preliminary analysis 
suggests that the large-scale cluster temper- 
ature structure (particularly for z > 0) is a 
potential discriminator of cosmological model. 

Our current adiabatic catalogs of clusters 
were intended to serve as a well-defined base- 
line sample for statistical and individual com- 
parison with both observed clusters and with 
future simulations that include additional 
physical effects. Remarkably, our initial sim- 
ulations of clusters with more realistic physics 
(radiative cooling and galaxy feedback) give 
rise to temperature profiles (Fig. 6) that agree 
well with our adiabatic profile and suggest 
that this truly is a universal cluster tempera- 
ture profile. 
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